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We have experimentally investigated quantum interference corrections to the conductivity of 
graphene nanoribbons at temperatures down to 20 mK studying both weak localization (WL) and 
universal conductance fluctuations (UCF). Since in individual nanoribbons at millikelvin temper- 
atures the UCFs strongly mask the weak localization feature we employ both gate averaging and 
ensemble averaging to suppress the UCFs. This allows us to extract the phase coherence length 
from both WL and UCF at all temperatures. Above 1 K, the phase coherence length is suppressed 
due to Nyquist scattering whereas at low temperatures we observe a saturation of the phase coher- 
ence length at a few hundred nanometers, which exceeds the ribbon width, but stays below values 
typically found in bulk graphene. To better describe the experiments at elevated temperatures, we 
extend the formula for ID weak localization in graphene, which was derived in the limit of strong 
intervalley scattering, to include all elastic scattering rates. 

PACS numbers: 



I. INTRODUCTION 

Phase coherent effects in graphene are determined by 
the combined action of several scattering mechanisms. In 
the past, extensive studies have been performed on those 
effects in bulk grapheneir— . Little attention, however, 
has been paid to phase coherent behavior in graphene 
nanoribbons (GNRs) where lateral confinement causes a 
crossover from 2D to ID behavior and additional scatter- 
ing is introduced at the edges of the ribbons. 

In the experiments of Morozov et al£ on bulk graphene 
strong suppression of weak localization was observed. A 
theoretical description of the phase coherent effects was 
given by McCann et al£, where elastic scattering mech- 
anisms (intra- and intervalley scattering) determine if 
weak localization (WL), weak antilocalization (WAL) or 
none of them is observed. If there is neither intravalley 
scattering nor intervalley scattering weak antilocal- 
ization is found, as expected for chiral quasiparticlcs 
associated with Berry phase tt— . Intravalley scattering 
tends to suppress the chiral nature of quasiparticles and, 
thus, destroys localization, whereas intervalley scattering 
tends to restore the weak localization effect^. In further 
experiments the phase coherent effects could be inter- 
preted by this theoretical description^—. Furthermore it 
was found that by changing the carrier density and/or 
the temperature, it was possible to alter the ratio of 
various scattering rates and observe a transition from 
WL to WAL as the chiral nature of the charge carriers 
was restored*^. 

In the case of graphene nanoribbons, however, due to 
scattering at the edges, intervalley scattering is predicted 
to be the most important mechanism leading to the 
observation of weak localization^. To our knowledge, 
up to now there are no extensive experimental studies 
reported on the analysis of weak localization in GNRs 
and the theoretical predictions still need to be verified 
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FIG. 1: (a) Scanning electron microscope image of a typical 
sample from this work. The length of the GNRs is 1 /im, the 
width 70 nm. Two palladium contacts are visible, (b) Array 
of GNRs with a GNR length of 1 fim in between two palla- 
dium contacts, (c) Sample C: Two-terminal resistance as a 
function of Ve, s for different temperatures at zero magnetic 
field, (d) Zoom-in of the GNR array, every GNR has a width 
of 40 nm and a spacing of 30 nm to the next, the zoom-in 
area of panel (b) is marked in white. 



experimentally. 

Another correction to the conductivity are universal 
conductance fluctuations (UCFs), which appear when 
the sample length does not greatly exceed the phase co- 
herence lengths. In our samples they are clearly visible. 
In graphene, these fluctuations are sensitive not only to 
the phase coherence length and the thermal length, but 
also to elastic scattering (intervalley scattering and in- 



travalley scattering). The conductance variance strongly 
depends on the exact types of elastic scattering present 
in the sample and is a factor, a, times larger compared 
to a usual metal because of valley degeneracy. If all scat- 
tering effects are negligible a = 4. For weak interval- 
ley scattering and either strong intravalley scattering or 
strong trigonal warping a = 2 and for strong intervalley 
scattering a = l 11 ' 12 . Experiments on graphene analyzed 
the universal conductance fluctuations by the correlation 
functio n 13 ! 14 and showed that those fluctuations can be 
used, for example, for thermometry^ 3 -. But up to now 
no studies on one dimensional graphene structures were 
performed. 

Both effects, namely the weak localization as well as 
the universal conductance fluctuations, allow us to ex- 
tract the phase coherence length L v in an independent 
way. Therefore, we performed experiments on graphene 
nanoribbons to study both effects. 



II. EXPERIMENTAL DETAILS 

Single layer graphene is deposited on a highly doped 
silicon wafer with a 300 nm thick Si02 layer by conven- 
tional exfoliation^. The flakes were imaged under an op- 
tical microscope and their position was detected with re- 
spect to predefined markers. The graphene nanoribbons 
(GNRs), as well as the arrays of GNRs, were fabricated 
by electron beam lithography and oxygen plasma reactive 
ion etching. The ribbon length was I /jm and the ribbon 
width W varies between 40 nm and 80 nm. For the trans- 
port measurements palladium contacts were attached to 
the GNRs, using standard electron beam lithography and 
thermal evaporation. Micrographs of typical samples are 
shown in Fig. [T] Electronic characterization and mag- 
netotransport measurements were done in two different 
cryostats with temperatures ranging from 1.7 K to 125 K 
and T = 20 mK to 900 mK, respectively, with magnetic 
fields up to B — 16 T. The measurements were done 
in two terminal geometry using standard lock-in tech- 
nique with frequencies of 13 or 17 Hz and an excitation 
current of 10 nA at Kelvin temperatures and 0.5 nA at 
millikelvin temperatures, respectively. To induce charge 
carriers in GNRs a gate voltage up to ±80 V was applied 
between the graphene and the Si wafer, see Fig. QJc) for 
typical backgatc measurements. Conductance measure- 
ments as a function of backgatc voltage, temperature and 
magnetic field were done on many different devices yield- 
ing consistent results. We show representative data for 
two individual graphene nanoribbons and two graphene 
nanoribbon arrays. For sample parameters and studied 
temperature range see Table HI 



III. RESULTS AND DISCUSSION 

Figure[2ta) shows the magnetotransport data collected 
for a 40 nm wide individual GNR at temperatures from 




FIG. 2: (Color online) The conductance G as a function of 
magnetic field B of a individual GNR shows quantum inter- 
ference phenomena for temperatures from T = 1.7 K to 48 K 
(a) and T = 50 mK to 900 mK (b). 



T = 1.7 K to 48 K. Weak localization is observed at 
low fields (\B\ <1.5 T) as well as universal conductance 
fluctuations, whose amplitude increases with decreasing 
temperature. 

For mK- temperatures large universal conductance 
fluctuations overlay the weak localization feature of the 
individual GNRs [Fig. [2jb)]. In order to still determine 
phase coherent properties different methods can be used: 
(i) Gate averaging: by adding up the measurement traces 
of different gate values one obtains an average conduc- 
tance which shows a clear conductance dip and which 
allows to fit weak localization, (ii) Ensemble averaging: 
Structuring an array of many graphene nanoribbons in 
parallel suppresses the UCFs and the phase coherence 
length can be obtained from fitting the weak localization 
feature. 

Furthermore we analyze the universal conductance fluc- 
tuations: the phase coherence length can be determined 
by calculating the autocorrelation function of the UCFs 
or by analyzing the amplitude of the UCFs. 
In the following sections we present all the different meth- 
ods mentioned above and finally we compare the results. 
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FIG. 3: (Color online) Sample A: Conductance G as a func- 
tion of magnetic field B at (a) T = 1.7 K and (b) T = 48 K. 
The weak localization feature is fitted by using Eq. [T] (orange) 
and Eq.[2](blue dashed line). Both fit formulas reproduce the 
data well at low temperatures (a), but at higher temperature 
Eq. [2] is more appropriate. 



Weak localization in individual graphene 
nanoribbons 



First we analyze and interpret the weak localization 
effects of single GNRs by fitting to theory&il. The re- 
sistivity correction 6p(B)/p 2 is given by the following 
formula valid in the limit of a very short intervalley scat- 
tering time^ and by fitting the magnetotransport data 
one can determine the phase coherence length L v . 
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with the diffusion coefficient D, the dephasing 
time t v = L 2 V D~ X and the magnetic relaxation time 
t b = 3h 2 /(DW 2 e 2 B 2 ). The above formula is valid if the 
magnetic length L rn = y h/eB is larger than the rib- 
bon width. For some samples, L m is on the order of W 
in the field range considered here. However, since the 
phase coherence length is extracted from the behavior 
around B = 0, the ID formula is still appropriate to de- 
termine L v . Also, we fitted the data with both the ID 
and 2D formula^ and found that the 2D formula was not 
able to describe the data well. Having a closer look at 
the individual fits, one recognizes that at low tempera- 
tures this simple fit formula (Eq. [TJ reproduces the weak 



localization feature well at low temperature [Fig. |3]Ja)] . 
However at higher temperatures the magnitude of the 
effect is overestimated [Fig. EJb)], because the phase co- 
herence length and the intervalley scattering length are 
of the same order. 

Therefore we generalized Eq. Q] to account for a fi- 
nite intervalley scattering time by including other rel- 
evant elastic scattering times, as done previously^ for 
two-dimensional graphene. Usually, the WL correction 
is described in terms of particle-particle correlation func- 
tions, so called Cooperons. In two-dimensional graphene 
8g is determined by the interplay of one pseudospin 
singlet (Cq) and three triplet (Cq,Cq,Cq) Cooperons, 
Sg oc —Cq+Cq+Cq+Cq and their corresponding re- 
laxation rates (cf. Ref.— ). For graphene nanoribbons the 
four Cooperons Cq , Cq , Cq and Cq need to be considered 
in a similar fashion. Therefore we have to include the con- 
tributions from one Cooperon Cq (with 2 , where Tj 
is the intervalley scattering rate) and from two Cooper- 
ons Cq and Cq (with r" 1 , which includes both the inter- 
and intravalley scattering rates). This leads to the fol- 
lowing formula^: 

-1/2 
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Here all scattering terms relevant in two-dimensional 
graphene (r" 1 , Tg , t~ 1 and r" 1 ) are included, with 
the corresponding lengths = \/Dt v ^^. Fitting 

the data with Eq. [2J with the intervalley scattering 
length Li about the ribbon width and the inter- and 
intravalley scattering length L» about a few nanometers, 
one obtains much better fits than with Eq. [T] especially 
for higher temperatures, cf. Fig. [3J For sample A a 
phase coherence length L v between 50 nm and 100 nm 
can be extracted, Fig. |6ja). As it turns out, the phase 
coherence lengths obtained by Eq. [T] and Eq. [3] are very 
similar, which proves the robustness of L v and confirms 
the validity of Eq. [2] 

Decreasing the temperature, universal conductance 
fluctuations strongly overlay the WL feature, cf. 
Fig. [U[b). In order to extract the weak localization and 
thus the phase coherence length one can do an averaging 
over different gate voltages. For sample B, the magneto- 
transport was measured for different gate voltages (from 
-40 V to -20 V in steps of 1 V) at T = 20 mK. The 
individual traces show strong conductance fluctuations, 
but by adding up those twenty-one measurements one 
obtains an average conductance which shows a clear con- 
ductance dip, Fig. SI Fitting this feature one obtains a 
phase coherence length of 100 nm [Fig.[l){c), purple star], 
which is in a reasonable order of magnitude. 
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FIG. 4: (Color online) Sample B, gate-averaging: The mag- 
netic field dependence of the conductance was measured at 
different gate-voltages and the arithmetic mean was calcu- 
lated. The average conductance G of the 40 nm GNR at 
T — 20 mK clearly shows the weak localization feature. Fit- 
ting the conductance dip by Eq. [T] and Eq. [2] yields a phase 
coherence length of 100 nm. 



We also performed a number of gate dependent 
experiments on some of the samples but we only saw 
a decrease of L v around the charge neutrality point, 
cf. RefJ^. We never observed a transition to weak 
antilocalization. Following Tikhonenko et al£ this is also 
not to be expected since they observed the transition in 
clean bulk graphene when t v became shorter that the 
clastic scattering times. Here, due to nanopattcrning, 
the elastic scattering times are so short that we could 
not reach this regime^. 
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FIG. 5: (Color online) Magnetotransport data of a GNR ar- 
ray. In comparison to the data of individual ribbons Fig. [2] 
the measurements of the arrays clearly show a suppression 
of the universal conductance fluctuations for all temperatures 
(a) and (b), whereas the weak localization feature is not af- 
fected. Blue dashed lines are best-fit curves to Eq. [2] 



C. Conductance fluctuations 



B. Weak localization in arrays of graphene 
nanoribbons 

Also ensemble averaging by measuring arrays of 
graphene nanoribbons suppresses the UCFs. Therefore 
arrays of graphene nanoribbons were fabricated and the 
conduction per ribbon was calculated. As expected, 
the parallel arrangement of the nanoribbons leads to a 
suppression of the universal conductance fluctuations, 
whereas weak localization is not suppressed [Fig. Eta) 
and (b)]. Thus the phase coherent effects can be 
separated and the weak localization feature can be fitted 
again with Eq. [1] or Eq. [2] which were introduced in 
section III A for individual graphene nanoribbons. 
Fitting the WL dips for sample C and D [Fig. [5Ja) and 
(b)] to Eq. Q] and Eq. [21 a phase coherence length L v 
between 30 nm and 80 nm can be extracted for the 
array of sample C and between 80 nm and 170 nm for 
sample D, Fig. Etb) and (d). 



Other methods to determine the phase coherence 
length base upon the analysis of the universal conduc- 
tance fluctuations. Therefore let us first interpret the 
data via the autocorrelation function^. The correlation 
of the conductance fluctuation is given by the correlation 
field Be and can be determined from the autocorrelation 
function defined by F G (AB) = j dB G(B) ■ G(B + AB), 
where the integration was done at magnetic field ranges 
not including the weak localization feature. The correla- 
tion function is normalized to the value at B = T and 
the correlation field is defined thus that the function 
drops to half the maximum value, Fq{B c ) = 0.5 • Pb(0). 
The phase coherence length can be extracted from 
the correlation field: L v = C\ ■ B c ■ W/&o, with C\ 
a prefactor between 0.95 for ^> Lt and 0.42 for 
L v <C Lj^ and <£>o the magnetic flux quantum. We 
determined the phase coherence length L v for different 
temperatures giving values between 100 nm and 500 nm 
for an individual ribbon (sample B) and between 130 nm 
and 530 nm for the array of GNRs (sample D). The 
values of L v for sample B and D are summarized in 
Fig. Oc) and (d), respectively. 
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FIG. 6: (Color) Phase coherence lengths determined via different methods for (a) sample A, (b) sample C, (c) sample B and 
(d) sample D. The phase coherence length L v was determined by different methods like weak localization (WL1 and WL2), the 
amplitude of the UCFs (RMS) and the autocorrelation function (AUT). The data points obtained by fitting the WL feature 
to the simple formula (Eq. [1]) are plotted as black, open squares (WL1) and by fitting the full formula (Eq. [2]) as black, filled 
squares (WL2), respectively. The orange, dashed lines represent the thermal length Lt- 



The rms amplitude AG rms of the fluctuations allows us 
to extract the phase coherence length in an independent 
way. From the WL measurements we conclude that the 
Ltp value is lower than the length of the ribbon L. The L v 
value can be less than or greater than the thermal length 
Lt = (Dh/ltBT) 1 / 2 , where D is the diffusion constant. If 
L v < L, Lt then AG rms depends on L v by the following 
relation for ID: 

e 2 (L \ 3/2 
AG rms = a-C 2 - l-f) , (3) 

where we set a = 1 due to strong intervalley scattering 
that mixes the valleys completely, as already shown by 
the analysis of the weak localization. In the temperature 
range considered here, the prefactor Ci ranges from 1.6 
to 2.4i£. In Fig. [SJc) the temperature dependence of L v 
extracted from the UCFs obtained by sweeping either 
the magnetic field (red) or the back gate voltage (blue) 
is shown for sample B. In contrast to Ref<22, we do not 
observe a breakdown of the ergodic hypothesis. Rather, 
as expected, we find similar fluctuation amplitudes 



of about 0.4 e 2 /h and thus the values of L v deduced 
from AG rms (Vb g ) and AG rms (B) match extremely well. 
Furthermore the absolute values of L v extracted from 
the rms amplitude AG rms and from the autocorrelation 
match very well and the temperature dependence is 

^ ji-0.19 

For GNR arrays of N ribbons the total conductance is 
given by Gjv = N ■ G\, with G\ the conductance of a 
single ribbon. The absolute conductance of the array 
is ~ N times larger than for an individual GNR. For 
further analysis, the variance of the conductance is 
calculated as var(GAr)= N- var(Gi) and the average 
fluctuation amplitude as AGjy = VN ■ AGi. Thus 
ensemble averaging increases the conductance amplitude 
only by a factor ~ \/~N and for the determination of L v 
by analyzing AG rms of graphene nanoribbon arrays one 
has to take this factor into account 2 ^. 
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D. Discussion 

Fitting weak localization with the standard fitting 
formula for narrow wires (Eq. [T]) was appropriate only 
at low temperature, but reaches its limit of applicability 
at Kelvin temperatures. Therefore we expanded the 
standard formula to Eq. [2] At high temperature, the 
corresponding fits describe the measured data much 
better. At mK-temperatures universal conductance 
fluctuations mask the WL feature. Different averaging 
methods (gate- and ensemble-averaging) allow us to 
still analyze the sample properties. Furthermore the 
amplitude and the autocorrelation function of the 
universal conductance fluctuations themselves were 
analyzed. The determined phase coherence lengths are 
comparable to the values of L v obtained by fitting the 
weak localization, whereas the values of L v determined 
by the autocorrelation function are always slightly higher 
than those obtained from other methods. We note that 
there is a discrepancy between the values of L v obtained 
from the autocorrelation and the rms amplitude of the 
UCFs for the arrays of nanoribbons [Fig. [6jd)] . This 
may be due to bulk graphcne leads, contributing more to 
the total resistance than for a single GNR, thus making 
the analysis of the correlation field less reliable. 

Figure [6] summarizes the values of the phase coherence 
lengths determined by different methods. Theoretically 
the phase coherence lengths is determined from the phase 
coherence time t v by L v = Dt v . For Nyquist scat- 
tering Tip is proportional to T -2 / 3 and therefore L v is 
expected to be proportional to T~ x / 3 ——. In our ex- 
periment, the temperature dependence of L v is about ~ 
T -0 - 3 at Kelvin temperatures and hence agrees with this 
model. But for mK-temperature it gets weaker, suggest- 
ing a saturating behavior at a few hundred nanometers. 
Thus L v clearly exceeds the ribbon width for most of our 
samples, suggesting that the etching process (at the GNR 



fabrication) does not severely reduce the phase coherent 
properties of the sample. However, the values of values of 
L v in graphcne nanoribbons (and graphene antidot lat- 
tices 2 ^) seem to be smaller than in bulk graphene. This 
could in principle be a consequence of the reduced diffu- 
sion constant D. For our GNR samples A to D we obtain 
T v,A= 41 ps, t^b= 3 ps, t V}C = 0.9 ps and t Vi d= 3.8 ps 
at the lowest temperatures while a bulk graphene sample 
with D = 0.046 m 2 /s showed L v = 2 /im and t v = 100 ps 
at 300 mK. We conclude that for all nanoribbon samples 
the phase coherence time as well as the phase coherent 
length are smaller than in bulk graphene. 

The presence of spin flip processes^ might explain the 
lower values of L v in graphene nanoribbons compared 
to bulk graphene 2 ^: Localized spins at the ribbon 
edges may lead to a de-phasing by spin flip scattering 
and thus lower the phase coherence length. Having 
localized spins at the ribbon edges, one could think of ex- 
periments with graphcne nanoribbons as spin injectors 2 ^. 

IV. CONCLUSIONS 

In conclusion, we have performed magneto-transport 
measurements in graphene nanoribbons as well as in 
arrays of GNRs. The observation and analysis of weak 
localization and universal conductance fluctuations allow 
us to determine the phase coherent properties of those 
graphene nanostructures. 
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Sample N W(nm) Temperature fj, (cm 2 /Vs) n (10 16 m~ 2 ) L mJp (nm) D (m 2 /s) 



A 


1 


40 


1.7 K - 


125 K 


165 


5.3 


4.4 


0.022 


A 


1 


40 


20 mK - 


900 mK 


165 


5.3 


4.4 


0.022 


B 


1 


40 


20 mK - 


900 mK 


330 


50.6 


27.4 


0.014 


C 


46 


40 


1.7 K - 


125 K 


680 


5.4 


18.5 


0.009 


D 


23 


80 


20 mK - 


900 mK 


500 


4.3 


12.1 


0.006 



TABLE I: Characteristic parameters of the different samples: Number of ribbons N, ribbon width W, mobility /i, carrier 
density n, mean-free path L m f P and diffusion constant D. 



